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O ' I. INTRODUCTION 

(N 

The recent progress in atomic clocks is largely due to improved methods for cooling, trapping and manipulating 
^ ^ ■ cold atoms. Improvements in the accuracy of atomic clocks continues unabated, and reports of progress in fountain 
clocks [1, 2] and optical lattice clocks [3, 4] are published monthly. An atomic clock with one part in 10 16 accuracy 
has been recently reported [5-7]; this is orders of magnitude better accuracy than the Cs standard atomic clocks. But 
collisional frequency shifts can significantly impact both the stability and accuracy of atomic clocks. These effects have 
i— - n , been systematically studied, e.g., for hydrogen [8], for Cs and Rb microwave clocks [10-12] as well as Sr optical clocks 
[5-7, 13]. It was conclusively demonstrated in the experiments reported in Refs. [6, 7] on a fermionic Sr atomic clock 
Oh | that quantum statistics can play a critical role in shaping interactions between atoms, and therefore in determining 
-t— > ■ collisional clock shifts. When the temperature of the fermionic atomic cloud is so low that the gas is degenerate, the 
Pauli exclusion principle forces atoms in the gas that are in the same internal state to occupy different vibrational 
levels, since their spatial wave function must be antisymmetric. These atoms cannot s-wave scatter and therefore, 
qh at low temperatures, where higher partial waves are frozen out, no collisional shift is possible. Refs. [6, 7] showed 
that a small component of the probe beam along the weak confining direction leads to different Rabi frequencies for 
the weak confining direction modes, and the optically induced inhomogeneity transfers some of the atoms into an 
| antisymmetric internal state, and this antisymmetric state can s-wave scatter. This gives rise to a small collisional 
shift. 

The analysis of the experimental results in Refs. [6, 7] was carried out assuming that the internal states of two 
(T) . colliding atoms in the gas are in a pure state [^(\ipi}\il>2) — |V'2)|'0i)): as m Eq. (13), where \tpi) and \tp2) can be 
superposition states of the ground and excited clock states, as in Eq. (10)]. Moreover, reference [14], which introduces 
a microscopic many-body approach to treat collisional shifts, Ref. [15], which reports exact calculations of collisional 
frequency shifts for several fermions or bosons, and Ref. [16], which also calculates collisional shifts due to interatomic 
interactions, all assume a pure internal state for colliding atoms in their analysis of collisional shifts. Here, we 
generalize the theoretical treatment of exchange symmetry considerations in collisional clock shifts to include mixed 
(density matrix) internal states for two colliding atoms, pab, int [see Eq. (2)]. That is, our treatment does not assume 
that the colliding atoms are in a pure internal state, but allows for mixed states to be treated. We shall also apply our 
analysis to bosonic-atom atomic clocks. We assume that the atomic gas is sufficiently cold that only s-wave collisions 
■ can occur, with higher partial waves frozen out because of the angular momentum barrier, but we do not assume 
that the gas is fully degenerate. If the atoms in the gas are fermionic, but are not in the quantum-degenerate regime 
(i.e., is not described by a degenerate fermi sea of fermionic atoms), a mixed state description is in general required. 
Similarly, bosonic atoms that are not in the quantum-degenerate regime, require a mixed-state description. We use 
the mixed state (density matrix) description of Refs. [17, 18] to develop a formalism to treat the correlation and 
exchange symmetry considerations necessary to properly describe clock shifts due to collisions of atoms in a mixed 
internal state. We are therefore able to determine how the fact that the gas in Refs. [6, 7] is not quantum-degenerate 
affects the calculation of the clock shifts. As we shall show, having a mixed state, rather than a pure state for the 
internal state of the colliding atoms, does not strongly affect the calculated results for the experimental conditions 
described in Refs. [6, 7]. 

The density matrix formalism we develop here could also allow for calculation of blackbody radiation effects for 
atomic clock transitions. With the experimental demonstration of suppression of collisional clock shifts due to strong- 
interaction many-body effects [19] (which we do not attempt to model here, since this requires determining the 
correlated many-body state of the atoms in a strongly interacting gas, and is out of the scope of this work), calculation 
of blackbody radiation shifts becomes even more important, since this can be the largest remain contribution to the 
error budget of the atomic clock. 
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The outline of the paper is as follows. In Sec. II we discuss the effects of exchange symmetry on the collisional 
frequency shift in atomic clocks, Sec. Ill applies the concepts developed in Sec. II to a fermionic atomic clock. In 
Sec. IV we calculate the collisional shift for a 87 Sr atomic clock with parameters similar to those used in Refs. [6, 7]. 
Section V develops expressions for the clock shift of a bosonic atomic clock, and a summary and conclusion are 
presented in Sec. VI. 



II. EXCHANGE SYMMETRY CONSIDERATIONS 



The density matrix characterizing systems of identical bosonic or fermionic atoms must be properly symmetrized. 
This can be accomplished by applying symmetrization or antisymmetrization operators. For a degenerate gas of N 
bosonic atoms, the state can be symmetrized by applying the symmetrization operator <S, whereas for fermionic atoms 
it must be antisymmetrized by applying the antisymmetrization operator A. The symmetrization operator for an N 
particle state is S = J2p P/N\ and the antisymmetrization operator is A = J2p $p P/Nl, where P is a permutation 
operator, P = ( n ^ "' ^ ), the sum is over all the N\ permutations of the atoms, and 5p = (— l) p [20]. For example, 

for two atoms (say, that are undergoing a collision in a gas), S = ^(1 + Pab) and A = 2(1 — Pab)- Thus, properly 
symmetrized two-particle boson and fermion density matrices are of the form 

Pab = 5(1 + Pab)pab\{1 + Pab) , (la) 



PAB Sym = 5(1 " Pab)pab\{1 Pab) ■ (lb) 

The density matrix pab for two colliding atoms has both an internal part, pAB,mt (for atoms in an atomic clock, 
the internal part can be restricted to involve only two atomic levels, the ground and excited levels participating in 
the optical transition of the clock), and an external part, pAB,cxt, involving the motional degrees of freedom of the 
atoms: 

PAB = PAB.cxt ® PAB.int ■ (2) 

If the atoms arc free, or if they are confined by a harmonic potential, the external part can be separated into a center 
of mass and a relative motion part, 

PAB.cxt = PAS, cm <& PAS, rei ■ (3) 

A symmetric density matrix for the internal degrees of freedom, pab, int = Pab' must be multiplied by a symmetric 
[antisymmetric] spatial density matrix pab, rei in the spatial degree of freedom tab — *A — r B for bosons [fermions], 
whereas an antisymmetric internal density matrix, pas, int = Pab> must be multiplied by an antisymmetric [symmetric] 
spatial density matrix Pab, rei for bosons [fermions], so that the full density matrix has the right exchange symmetry. 
For cold gas temperatures, all but I = partial waves are frozen out due to the angular momentum barrier potential 
Vi(r) = h 2 l(l + l)/(2pr 2 ), where r = \?ab\ is the relative distance between the colliding atoms and p is their reduced 
mass. Since s-wave scattering wave functions are symmetric in the exchange of the colliding particles, we conclude that 
at cold temperatures, the internal density matrix must be antisymmetric for colliding fermionic atoms and symmetric 
for bosonic atoms. 

For the internal degrees of freedom of the colliding atoms, the symmetrization operator is S = ^(1 + Pab) = 
| + \cr a ■ <J b and the antisymmetrization operator is A = ^(1 — Pab) = t(1 — ""A ■ o"b) ; hence, the density matrix 
for the internal degrees of freedom are: 

Pab = hj + \°A ■ cr B j PAB f| + ^<ta • ersj , (4a) 



Pab = \ (1 - o A ■ <t b ) pab \ (1 - a A ■ tr B ) . (4b) 

Quite generally, the density matrix for the internal degrees of freedom of two two-level systems can be written in 
terms of the Pauli matrices, <x, for each of the two-level systems as follows [17, 18]: 

Pab = | [(1 + n A ■ cr A ) (1 + n B ■ er B ) + a A • C A b • vb] ■ (5) 
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Here rii = (er^) = Trp<Tj is the expectation value of the 'spin' (recall that this is the spin analogy used to describe 
our two-level system), as can be verified by multiplying (5) by er^ and tracing over the internal states of atoms A and 
B, and the correlation matrix is Cij_AB = (&i.ACj,B) — (ci,A)(o~i.B), as can be verified by multiplying (5) by ai^A^j.B 
and tracing. If the two colliding two-level atoms are in an uncorrelated (and therefore unentangled) state, Cab = 0, 
and the density matrix pAB reduces to a tensor product of the density matrices for each of the internal states of the 
atoms, pab = PA® PB- 



III. FERMIONIC ATOM CLOCK 



The trap potential holding the atoms is taken to be a harmonic oscillator potential (see Sec. IV), but at short 
range, the potential is the molecular potential of the atoms. Hence, the potential between two atoms in the trap is 
the molecular potential at short range (nm scale) and at large distance (/im or mm scale) it becomes the asymmetric 
harmonic potential confining the atoms. The external (trap) states of the atoms in the atomic clock will be taken to be 
thermally populated harmonic oscillator states at large relative intcrnuclcar distances (see Sec. IV), but in this section 
we concentrate on the internal state of the atoms. The internal state for fcrmionic atoms must be anti-symmetrized at 
low temperatures; hence, that part of an arbitrary initial density matrix p for two atoms that s-wave scatters is given 
by ApabA = >^Psingieti where /9 s i ng i t = j(l — er a ' ffl) is the density matrix of the singlet state and the correlation 
coefficient x = TyApabA = \ [1 - (n A ■ n B + J2i Cu)], i- e -> 



1 



% Psinglct 



1 - (nx- n B + ^ Cg) 



-(1 - a a ■ a- B ) 



(6) 



The quantity x = j [1 — (n^ -115 + Y^,i Cu)] is the fraction of the density matrix that is in the singlet state. Only 
this part of the initial density matrix is able to s-wave scatter when the atoms are fermions. The expression in Eq. (6) 
plainly involves the correlation of the spins, both because of the the spin-correlation term involving the C matrix and 
the riA ■ lis term. The clock shift for fermionic atoms is given by the expression [6, 8, 10, 12, 23, 24] 



1 
T 



dtAv(t) 



2h 



L s,ge 



T 



dtx(t) [nJt) - n e (t)} 
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where a s ge is the s-wave scattering length for collisions of atoms in states g and e, and the correlation coefficient 
x (t) is the instantaneous collisional shift parameter for fermions, equal to the the pair correlation function G 1 ^ 1 (0) 
for collisions of pure internal states, but more generally is given by 



x(t) 



1- (nA(t)-n B (t) + Y^Ca(t) 



(8) 



n g (t) and n e (t) are the densities of atoms in the ground and excited state respectively [n g = gp gg , n e = gp ee ] 
and g = N/V is the density of atoms in the gas. Note that the instantaneous collisional shift A^(t) = 
(2h/m)a St g e x(t) \n g {t) — n e (t)}, i.e., is proportional to x(t) [n g (t) — n e (t)]. 

If the two colliding atoms are in a pure internal state with the atoms in the same superposition state, IV'iIa = |V , i)b> 
then the two-particle internal state is a product state with • n B = 1 and Cab = 0, so x(t) = 0, and no collisional 
shift is present. Since, as explained in Refs. [6, 7], atoms in different trap states Rabi-flop at different rates, collisional 
shifts in pure-internal states are nonvanishing. It is instructive to consider the case of a two-atom density matrix that 
corresponds to a product state, 



I*) = \lk) A \lh)B , 

with if) 1 ^ ip2- This is an uncorrelated (unentangled) pure internal state with 

\^) A = a\g) + p\e) = (A , \^)b = j\g) + S\e) = ^ 



(9) 



(10) 



whose density matrix in (5) is also a product density density matrix, pab — PA® Pb, and Cab — 0. Hence, Eq. (8) 
becomes 



x(t) 



1 



(l-n A (t)-n B (t)) 



(11) 
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where the unit vector for particle j (— A, B) is = (sin 9 3 ; cos ■, sin 9j sin <j>j , cos 6j ) = (xj,yj,Zj). Up to a 
multiplicative phase factor, 

and similarly for \ip2)B — A simple example may be helpful. Suppose particle A is in state \g) = \\) and particle 
B is in state |e) = \\), and the particles are uncorrclatcd so pab — Pa® Pb', then = — z and tub = z, and using 
(11), x — 1/2. This result can be easily understood by noting that the state can be written as a fifty-fifty incoherent 
superposition of a singlet and triplet density matrix, hence x = 1/2. 

Let us compare with the case when the gas is fully coherent and the internal state of two fermions that interact via 
s-wave scattering is in an antisymmetric pure entangled state, 

\* A8 ) = ^(liMhM - KMVi}) • (13) 

This is the form of the internal state taken in the collisional shift calculation of Ref. [6]; the internal state of the 
two colliding particles is a fully coherent antisymmetric state (13), and therefore the internal state density matrix (5) 
corresponding to this state has only one non-zero eigenvalue, which is unity. The collisional clock shift Av is then 
proportional to 

k(t) = G^(0) = (#as\*as) = 1 - Hth*(t) +P(t)S*(t)\ 2 . (14) 

(2) 

x(t) for the pure state fermion case reduces to the general result x(t) in Eq. (8), i.e., for a pure state, x = x = G AS (0), 

(2) 

where G AS (0) is the pair correlation function for fermions at zero relative distance. For the specific simple example 
suggested above, a = 1, y3 = 0, 7 = and 6 = 1, so x = 1, whereas x = 1/2. Hence, the collisional clock shift for the 
mixed state is half that obtained for the antisymmetric state. Reference [6] assumes the internal states of the atoms 
in the gas is fully coherent, as in a fully degenerate gas, and therefore the internal quantum state of the two particles 
colliding in an s-wave collision must be antisymmetrized as in Eq. (13). We do not necessarily assume that the gas 
is degenerate, hence we use a density matrix description. In our simple example, the two colliding particles in the 
gas are either (with probability 50%) in an internal antisymmetric state (in which case the motional state must be 
symmetric) or in an internal symmetric state (in which case the motional state must be antisymmetric so only odd 
partial waves occur, but since the gas is cold, these collisions are frozen out, i.e., only s-wave collisions occur). Hence, 
in the simple example considered here, both possibilities occur with equal probability, and the factor x = 1/2. Note 
that, if the gas really is fully degenerate, only the antisymmetric internal state (i.e., the singlet) is populated, and x 
is be unity, as obtained in Rcf. [6]. 

The collisional shift parameter x for the state in the specific simple example above is the largest possible for any 
product state, and it has x = 1/2. The largest possible collisional shift parameter for an antisymmetric pure state is 
x = 1. The smallest collisional shift is for the product state with • = 1, i.e., \ipi) = \ip2)> f° r which x = 0. The 
smallest possible collisional shift parameter for an antisymmetric pure state is x — 0. For any given |?/>i) and \ip2), 
x for the product state is always smaller than k for the antisymmetric state (13), x < x. This suggests that it may 
not be worth going to the degenerate gas limit for a fcrmionic gas, i.e., it would suggest that the collisional shift is 
smaller if the gas is not fully quantum degenerate. 

Note that the uncorrelated case, pab = PA® Pb in (9), and the pure entangled state case corresponding to as) 
in Eq. (13), are the two extreme cases of states that are possible; mixed correlated and mixed entangled states fall 
in between these two limits [17]. Note also that our formulation allows us to specify any many-body state of the gas 
and thereby calculate the clock shift for arbitrary many-body gas states by tracing out all but two particles, thereby 
obtaining a two-particle density matrix. 



IV. MODELING THE 8 SR OPTICAL CLOCK 



In the experiments reported in Refs. [6, 7], a gas of 87 Sr 1 S$ {F = 9/2) atoms in a magic wavelength optical 
lattice [4] with wavelength A = 813.43 ran, which to a good approximation forms a harmonic trap with frequencies 
uj x = ujy = 2tt x 450 Hz and uj z «a 2tt x 80 KHz, is optically pumped to the tuf = 9/2 (or the mp = —9/2) state. 
The excited clock state, |e)(= |t)), is the 3 Po state and the clock transition has a wavelength of ~ 698 nm. The gas is 
at a sufficiently high temperature (e.g., T = 1 pK or 3 pK) that the harmonic oscillator motional states (n x , n y , n z ) 
of the atoms in the gas can be considered populated according to a Maxwcll-Boltzmann distribution (at T = 1 pK, 
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= n^ = 46). The laser probing the clock transition propagates along the ID lattice axis, and the atoms in the gas 
have a slightly different Rabi frequency Q nm ,n y ,n- depending on which trap state they populate. With perfect alignment 
of the probe laser along the strong confinement axis, a residual angular spread AO between the probe laser beam and 
the optical lattice k remains due to the finite size of the lattice beam. Moreover, if the symmetry is broken due to either 
aberrations in the beam profile or angular misalignment A9 between the lattice and the probe beam, an even larger 
effective misalignment results. As explained in Ref. [7], since the transverse trap is isotropic, a small net misalignment 
angle AO along x results in a wavevector of the probe laser that is given by k rj k(z + AO x). Hence, there are slightly 

different Rabi frequencies for the different trap states, ft nx ,n y .n, — e'~' ,I ~' ,1 '~'' i " 2 I I i I (»Ji)-£n t (1)2), where 
f2o is the bare Rabi frequency, L Ui are Laguerre polynomials, and rji are the Lamb-Dicke parameters for the transverse 

and longitudinal directions, Vx = r ]y = S^^^^II Vz = TZ^^l I 6 > 21 1- 

We now calculate the clock shift for such a cold but non-degenerate gas with a thermal population of trap states 
and an internal state so that two colliding atoms are in a non-correlated product internal state. The internal state 
density matrix is given at any time t by a general two-qubit classically-correlated state that takes the form [17] 



where k = n x ,n y ,n z , and P k = e~ Ek ^ T / ^2 k , e~ E *' l T is the Boltzmann distribution probability for state k (and 
similarly for Pi). Presumably {Pk} are constant in time. We used a Bolzmann distribution to calculate the probabilities 
for finding atoms in the various trap states in the calculations for the 87 Sr clock in Sec. IV A, since the experimental 
temperatures used were sufficiently high for this to be an excellent approximation. At lower temperatures we would 
use the correct fermi (or, in the case of bosons, bose) distribution. The evolution of the unit vectors iiA,k(t) and 
n B,k(t), and the population densities n g (t),n e (t) (which are also functions of the trap state k) are given by the 
evolution operator for the two-level system in the presence of the probe laser, U(t, 0) = e _lffLMt / s , where -Hlm is the 
light-matter Hamiltonian, is given by 



W^) + ffsin(^) f sin(^) 



We defined the generalized Rabi frequency, ft g = vT^wwna I 2 + an d the laser detuning A = lu — (E e — E g )/K. 
For an atom that at t = is in the ground electronic state and in trap state n x , n y , n z , the amplitudes a(t) and j3{t) 
of (12) are 

, , ,Q a t, *A . ,Q, a t. _ , . iQ, . Mat. , „, 

an.,n sl n.W = COs(^-) - — Sin(^-), P n *, ny ,nM = f^ Sln ( 2 ^ ( ^ 

The density matrix for the classically-correlated state can be written in the form of Eq. (5) with Bloch vectors 

n A(t) = Pk nA - k (*) ' nB (*) = H Pi ' ( 18 ) 

k I 

and correlation matrix vanishes, Cy = 0. The average clock-shift coefficient >?(£), where the average is over the trap 
states, is given by 

m = j2PkP l 1 - nAMt ^ nBAt) = \\i -j2wcose kA t)} , d9) 

k,l k,l 

where Qk.i(t) = riA,k(t) ■ riB.i(t). Note that H{t) cannot be written as a product of averages over the first and second 
particles separately since ®k,l(t) generally does not factorize into a product of terms depending separately on the first 
and second particle. In the clock shift calculation we need to average over motional states as follows: 

1 r T 2ft 1 f T 

Av = - / dt Av{t) = —a s . ge - / dt x(i) [n g (t) - n e {t)} , (20) 
1 Jo m 1 Jo 

where 

x(t) [n g (t) - n e {t)} = X - ^ P k [n 9tk (t) - n e , k {t)} j^P* [1 - n A , fc (i) • n B ,i(t)] | ■ (21) 
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Note that the quantity x{t) [n g (t) — n e (t)] is proportional to the average instantaneous frequency shift Av(t), i.e., 
x(t) [n g (t) - n e {t)} = ^L-A^(t). 

We also calculate the collisional clock shift for a gas with thermal population of trap states and an internal state 
so that two colliding atoms are in an entangled internal state with collisional shift parameter k(t), 

sw = E PfeP 'E PfcF '[ 1 -i a ( t )^w+/ 3 w r wi 2 ] • ( 22 ) 

k,l k,l 



The instantaneous collisional clock shift in this case is proportional to k[t) [n g (i) — n e (i)]. We are therefore able to 
compare and contrast this with Jc(t). Note that the average kit), can be written as an average of k(t) over n x , n y , 
n z for each of the two particles A and B, but this average also cannot be expressed by taking averages over a, f), 
and 7 and S separately because k(t) does not factorize into a product of terms depending separately on the first and 
second particle. 



A. Results of Sr Atomic Clock Calculations 



We now present numerical results that highlight the difference between using the two extreme limits for describing 
the internal states of two colliding atoms in the gas: a product state given by pab = Pa®Pb and a pure antisymmetric 
state as given in Eq. (13). 

In our calculations we use parameters for the system based upon those in Refs. [6, 7]. In all of our calculations we 
use a Rabi frequency Qq = 2n x 59 Hz. Figure 1 shows the average population inversion [n g (t) — n e (t)] versus time 
for detunings A = 2irx 0, 25 and 50 Hz for a gas at temperature T = 1 pK and for a laser misalignment of A8 = 10 
mrad. The atoms start off at t = in the ground electronic state and Rabi flop with time. For zero detuning, the 
oscillation of the average population inversion with time is about zero population inversion, but for finite A the time 
averaged population is different from zero, and grows with detuning since less atoms are put into the excited state for 
higher detuning. The decay of the Rabi oscillations with time results because of the thermal distribution of the atoms 
and the fact that the Rabi frequencies ^n x ,n y ,n z depend on the trap state populated in the thermal distribution. 




I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 

10 20 30 40 50 60 70 80 
time (ms) 



FIG. 1: Average population inversion, [n g (t) — n e (t)], versus time for three different detunings for T = 1 pK. The fact that the 
Rabi frequency depends on the motional state causes the amplitude of the Rabi oscillations to decay with time. 



Figure 2 compares the average correlation coefficient Tc{€) (heavy curves) and the pure-state k(t) (dashed curves) 
versus time for A = 27rx 0, 25 and 50 Hz and for a gas temperature T = 1 pK. The largest difference between H(t) 
and k(t) occurs for A = 0, and the difference diminishes with increasing A. The average correlation impacts the 
average collisional frequency shift [see Eqs. (19), (20) and (21)]. We see from the figure that these quantities increase 
with the clock run-time. For small detunings, the average pair correlation function at zero particle separation for 
pure states, kit) = G( 2 )(0), is substantially larger than x(i) for the product states, and the difference between these 
quantities diminishes with increasing A. This, of course affects, the collisional clock shift, as we shall see in the next 
figure. 

Figure 3 shows the quantity x{t) \n g {t) — n e (t)] versus time calculated for three detunings, A = 27rx 0, 25 and 
50 Hz, and for T = 1 pK. This quantity is proportional to the instantaneous clock shift Av(t) (but since the s-wave 
scattering length a Sjffe is not well known, the factor (2h/m)a Stge in Eq. (7) is removed from the quantity plotted). It 
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10 20 30 40 50 60 70 80 



time (ms) 

FIG. 2: (color online) 5r(t) (heavy curves) and the pure-state Zc(t) (light curves) versus time for three different detunings and 
for T = 1 /jK. 



is clear from Fig. 3, and from Eq. (7), that the temporal average of the clock shift for A = will be close to zero if 
the temporal average is taken over a time tf = T, which corresponds to a full Rabi cycle, but for finite detuning, the 

temporally averaged clock shift grows with detuning. 

Figure 4 compares the product state x(t) [n g (t) — n e (t)] (solid red curves) and the pure-state k(t) [n g (t) — n e (t)) 
(dashed black curves) for A = and A = 2tt x 50 Hz. The ps 25% difference between x(t) and k{t) evident in Fig. 2 
for A = at large times shows up mostly at the extrema of x(t) [n g (t) — n e (t)] and k(t) [n g (t) — n e (t)] in Fig .4; 
elsewhere the population-difference [n g (t) — n e (t)] becomes small and the difference between the averaged quantities 
plotted in Fig. 4 become small and h ard to s ec. As the detuning increases, the magni tude of t he co llisional phase 
shift is greater for k(t) [n g (i) — n e (t)] than for x(t) \n g (t) — n e (t)], but the differences in x(t) and k(t) decrease with 
increased detuning, as evident from Fig. 2. Thus, the non-degenerate gas has a lower clock shift (in magnitude) than 
the gas degenerate with an antisymmetric internal state. The effect appears to be very small, but in fact, we shall 
see in Fig. 8 that, partly due strong cancellation of the positive and negative clock shifts [particularly for A = - see 
Fig. 4(a)], the time averaged clock shifts x(tf) [n g (tf) — n e (tf)] and k(tf) [n g (tf) — n e (tf)] are not inconsequential. 




I 1 1 1 1 I 

10 20 30 40 50 60 70 80 
time (ms) 



FIG. 3: (color online) x(t) [n g (t) — n e (t)], which is proportional to the clock shift Au(t), versus time for three different detunings 
and for T = 1 fiK. 



Figure 5 compares x(t) [n g (t) — n e (t)] (solid red curves) and k(tf) [n g (t) — n e (t)] as a function of time for three 
different detunin gs at a relatively h igh temp erature of T = 10 / iK. The magnitude of the collisional phase shift is 
a little larger for k(t) [n g (t) — n e (t)] than for x(t) [n g (t) — n e (t)]. The difference increases with increasing detuning, 
but then saturates. The behavior here for T = 10 is actually quite a bit different from that shown in Fig. 4 
for T = 1 /xK. For example, the decay of the clock shift oscillations as a function of time at A = in Fig. 5 is not 
evident in Fig. 4, and the oscillations with time at A = 27rx50 Hz in Fig. 4 are largely damped in Fig. 5. Here (and 
in Fig. 4), for zero detuning, the instantaneous clock shift oscillates around zero with time but this changes with 
increased detuning. Note that the ordinate scales of the three fra mes in Fig. 5 arc different. 

Figure 6(a) plots the average excited state population n e (tf) and Fig. 6(b) the average clock shift 
k(t) [n g (tf) — n e (tf)] versus detuning A for tf = 80 ms (note that for Figs. 6-8, the bar indicating average also 
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FIG. 4: (color online) (a) >c(t) \n g (t) — n e (t)] (solid red curves), which is proportional to the clock shift Av(t), and the pure- 
state k(t) [n g (t) — n e (t)] (dashed black curves) versus time for detuning A = 0, and (b) A = 2ir x 50 Hz. The gas temperature 
is taken to be T = 1 fiK. 

means average over the clock time tf). The population for A = and tf = 80 ms is about 60% in the excited 
state, so n e (tf) > n g (tf). This is why the quantity x(tf) [n g {tf) — n e (tf)] is negative for A = 0. At the detuning 
A ?=! 2ir x 0.15 Hz for which n e (tf) = 0.5, clearly, k(tf) [n g (tf) — n e (tf)] = 0. At higher detunings, the oscillations in 
>c(tf) [n g (tf) — n e (tf)] as a function of detuning is due to the oscillations in the population n e (tf) versus A. 

Figure 7 plots n e (tf) and x(tf) [n g (tf) — n e (tf)} versus detuning A for T = 1 /iK and T = 3 /iK at a pulse duration 
of tf = 1.7 ms, which is near the start of the first Rabi cycle — see Fig. 1. Both the average excited state fraction 
in Fig. 7(a) and the collisional shift in Fig. 7(b) are smaller for the larger temperature. The collisional shift is seen 
to be very substantially smaller at T = 1 /xK than for T = 3 fj,K in Fig. 7(b). The temperature effect is dramatic for 
all detunings, i.e., k{tf) [n g (tf) — n e (tf)] is significantly reduced as the temperature is decreased from T = 3 /iK to 
T = 1 /iK for virtually all detunings. At A = 0, x{tj) [n g (tf) — n e (tf)] is positive [in contrast with Fig. 6 for tf = 80 
ms], since this is at the start of the first Rabi cycle (see Fig. 1) and n e (tf) <^ n g (tf). Moreover, very little vestige 
remains of the oscillat ions with detuning that was so prominent in Fig. 7. It is interesting to note that if one plots 
>c(tf) [n g (tf) — n e (tf)] versus n e (tf) [i.e., if one plots the ordinate of Fig. 7(b) versus the ordinate of Fig. 7(a)], the 
collisional shift is nearly linear with n e (tf) for both T = 1 /iK and T = 3/xK, but the slope is strongly temperature 
dependent and increases with increasing temperature. Holding n e (tf) (or the excitation fraction) constant, and 
lowering the temperature of the gas, serves to significantly lower clock shifts in this regime of low excitation fraction. 

Figure 8 plots the clock shift Av for tf = 80 ms as a function of temperature for four different detunings. For A = 
and A = 2ir x 15 Hz, the collisional shift first decreases with temperature from zero to around a value of —0.05, and 
at around 0.75 /iK, Av begins to increase with temperature. At large temperatures, Av appears to tend to zero. For 
the higher detunings of A = 2-k x 25 and A = 2-k x 50 Hz shown in the second panel, the c ollisional shift increases 
with increasing temperature and then saturates. Again, the magnitude of x(t) [n g (t) — n e (t)] is almost always smaller 
than k(tf) [n g (tf) — n e (tf)}. Moreover, we see a very different temperature dependence of the clock shift for small 
detunings than for larger detunings. 
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FIG. 5: (color online) x(t) [n g (t) — n e (t)] (solid red curves) and x(t) [n g (t) — n e (t)\ (dashed black curves) versus time for 3 
detunings, A = 0, A = 2n x 25 and A = 2n x 50 Hz, and for T = 10 ^iK. Note the different ordinate scales. 



BOSONIC ATOM CLOCK 



Let us now consider a non-degenerate gas of bosonic atoms that is nevertheless sufficiently cold so that only s- 
wave collisions participate in the collisional frequency shift. For example, consider an atomic clock based on neutral 
bosonic atoms trapped in a deep "magic-wavelength" optical lattice, or more specifically, 88 Sr atoms in the 5s 2 1 So 
\F = 0,Mf = 0). A weak transition to the 3 Po \F = 0,Mp = 0) state is available when a static magnetic field 
enables a direct optical excitation of forbidden electric-dipole character, that is otherwise prohibitively weak (the 
transition occurs by mixing the 3 Pi with the 3 Po state) [22]. Such an atomic clock, probed using Ramsey fringe 
spectroscopy [26], was analyzed in Ref. [27]. For the collision of two bosons in such a gas via s-wave collisions, the 
internal two-particle state must also be symmetric. The internal states of the atoms in such a collision can be of the 
form WD + \ID or ITT) i where \D = \g) and W) = \ e )- The s-wavc scattering lengths for these states arc denoted 
go and a Stee respectively. The Hamiltonian is given by [27], 



a\a e 



a\a g ) 



a\a g ) + 



Airti 2 



i,j=g,e 



(23) 



This leads to the following formula for the collisional frequency shift for bosons [12, 23, 24]: 



2% 1 f T ( 

Av = — — / dt I a s , oe ^ it4->+n.t) (t) \n (t) - n e (t)] + a s , gg K\u) (t)n a (t) + o a , ee J%f) (t)n e (t) } . (24) 
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This expression can be compared with Eq. (7) for the fermion case. Here, clock shifts are due to gg s-wave collisions, 
ee s-wave collisions and ge s-wave collisions. The three h correlation parameters that multiply these scattering lengths 
can be determined by calculating the projection of the density matrix of two colliding particles onto the symmetric 
states, ^= (|t|) + lit)'), Itt) and 14-1) ■ For example, 

TV) PAB(t) P(tt> = *Ttt) (t)P\tf) (*) i ( 25 ) 

where V^) = [ft") (ft I is the projector onto state |tt), and 

>c n) {t) = Tr^) pAB{t)V^)) . (26) 

After some algebra, we find that the correlation coefficients arc given in terms of the parameters of the density matrix 
in Eq. (5) by 

x tm+M) — T [1 + n A.x*iB,x + n A ,yn B .y - nA, 2 n B!2 + C xx + C yy - C zz ] , (27) 

V2 4 

*%\) = T [1 - n A,z ~ n B . z + n AyZ n B . z + C zz ] , (28) 
^|tf) = ^ [1 + "a, z + n B , 2 + n^^ns,* + . (29) 

When the atoms are not correlated, the Cy terms vanish in these formulas, and the correlation coefficients are given 
slowly in terms of the expectation values {nj,i} of the "spins". 

In order to use these results to calculate the clock shifts for atomic clocks using such atoms as 133 Cs and Rb 
and 88 Sr, three s-wave scattering lengths, a s gg , a s ^ ee and a s ,ge would are required. Hopefully, improved experimental 
values for these s-wave scattering lengths, perhaps in conjunction with theoretical calculations, will allow comparison 
of theory with collisional clock shifts for bosonic atomic clocks. 



VI. SUMMARY AND CONCLUSION 



In summary, we developed a formalism that can be used to calculate atomic clock shifts for arbitrary mixed states 
of an atomic gas. We applied the formalism to the fermionic 87 Sr atomic clock reported in Rcfs. [6, 7]. When the 
two fcrmions that collide in the gas are in a pure antisymmetric internal state, \^as) = ^(l^i)!^) — | ^2) |"0x ) ) , the 
formalism reduces to the one presented in Ref. [6]. For sufficiently high temperatures (e.g., T = 1 or 3 ^K), the gas 
may be closer to an uncorrelated gas, and we can treat this case, as well as the case of an arbitrary initial correlated 
state by projecting out the two-body density matrix and using it to describe the collision process of the atoms in the 
atomic clock. We saw in Sec. IV A that, for the experimental conditions described in Refs. [6, 7], when the internal 
state of colliding pairs of atoms in the gas is taken as product mixed state rather than an antisymmetric pure state, 
the collisional cloch shift is somewhat reduced, but the reduction is small. Moreover, we showed how the formalism 
can be used to describe a bose gas; then, all the s-wave scattering lengths ( 

^s,gg7 ^s,ee? Q*s, eg) are required to calculate 

the collisional clock shift [27]. 

The formalism developed here can also be applied to atomic clocks that use the Ramsey fringe double resonance 
technique [26, 27]. Moreover it can be used for clocks that have a strongly correlated initial state, since arbitrary 
two-particle density matrices can be modeled using the form in Eq. (5). Furthermore, even if the initial gas is a highly 
correlated many-body state, the two-body density matrix can be projected out of the many-body density matrix. The 
formalism might also be of utility for treating magnetometry wherein the collisions affect the spin-relaxation time 
depolarization via depolarization. 

Another possible application of the density matrix methods introduced here is to calculation blackbody radiation 
(BBR) shifts. BBR shifts in atomic clocks are caused by perturbation of the atomic energy levels of the ground and 
excited clock levels by the oscillating thermal radiation. Since we can carry out calculations with arbitrary density 
matrix (mixed) states, our method could be used to calculate blackbody radiation (BBR) effects, which can populate 
both ground and excited clock states of an atomic clock operating at microwave frequencies, if the clock operates 
at room temperature. Moreover, both ground and excited atomic levels of an optical frequency clock transition are 
perturbed by BBR, e.g., for the 1 5'o — > 3 Po clock transition of strontium, both clock levels experience a BBR shift 
[25]. The overall BBR correction is the difference of the BBR shifts for the two levels. If there is a small residual 
magnetic field present, the BBR could shift the energy of ground hyperfine levels due to coupling between the various 
F, Mp levels of the ground state. Such shifts would be in addition to shifts due to coupling of the ground and excited 
clock states to other states. A density matrix method of the type introduced here could be used to deal with the 
former type of clock shift. 
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FIG. 6: (a) Average excited state population n e {tf) versus detuning A for tf = 80 ms. (b) x{tf) \n g (tj) — n e (tf)] (solid red 
curve) and Zc(tf) [n g (tf) — n e (i/)] (dashed black curve) versus detuning A for tf = 80 ms. 




FIG. 7: (color online) (a) n e (i/) and (b) >c(tf) \n g (tj) — n e (tf)] versus detuning A for tf — 1.7 ms. 
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FIG. 8: (color online) *c(tf) \n g {tf) — n e (t/)] (solid curves) and Zc(tf) [n g (tf) — n e (tf)] (dashed curves) versus temperature for 
tf — SO ms and for four 4 detunings, (a) A = and A = 2ir x 15 Hz, and (b) A = 2ty x 25 and A = 2ir x 50 Hz. 
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